QCD dynamics in mesons at soft and hard scales 
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Abstract. Using a ladder-rainbow kernel previously established for the soft scale of light quark 
hadrons, we explore, within a Dyson-Schwinger approach, phenomena that mix soft and hard scales 
of QCD. The difference between vector and axial vector current correlators is examined to estimate 
the four quark chiral condensate and the leading distance scale for the onset of non-perturbative 
phenomena in QCD. The valence quark distributions, in the pion and kaon, defined in deep inelastic 
^7* scattering, and measured in the Drell Yan process, are investigated with the same ladder-rainbow 

^5 truncation of the Dyson-Schwinger and Bethe-Salpeter equations. 
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A great deal of progress in the QCD modeling of hadron physics has been achieved 
through the use of the ladder-rainbow truncation of the Dyson-Schwinger equations 
(DSEs). The DSEs are the equations of motion of a quantum field theory. They form 
^ an infinite hierarchy of coupled integral equations for the Green's functions (n-point 

'""' functions) of the theory. Bound states (mesons, baryons) appear as poles in the appro- 

, — I priate Green's functions, and, e.g., the Bethe-Salpeter bound state equation appears after 

r>" taking residues in the DSE for the appropriate color singlet vertex. For recent reviews 

— 1 on the DSEs and their use in hadron physics, see Refs. [1, 2, 3, 4]. 

cn In the Euclidean metric that we use throughout, the DSE for the dressed quark 

^ propagator is 

O /-A X' 

O S{py^ = Z2/i? + Z4m(Ai)+Zi / g^D^y{p-q)-Y^S{q)r^{q,p), (1) 
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where D^v(^) is the renormalized dressed-gluon propagator, T\{(i,p) is the renormal- 



;_( ized dressed quark-gluon vertex. We use / to denote / <i^/(2;r) with A being the 

"^ mass scale for translationally invariant regularization. The renormalization condition is 

^{p)~^ =17- P + tn{ii) at a sufficiently large spacelike /i^, with m{ii) the renormalized 
mass at renormalization scale /i. We use /i = 19GeV. The Zi{ii,A) are renormaliza- 
tion constants. Bound state pole residues of the inhomogeneous Bethe-Salpeter equation 
(BSE) for the relevant vertex, yield the homogeneous BSE bound state equation 



r^Hp+,P-) = rK{p,q-P)S^{q+)Y^^\q+,q^)S\q^) , (2) 
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where K is the renormalized qq scattering kernel that is irreducible with respect to a pair 
of qq lines. Quark momenta are ^+ = ^ + r\P and q^ = q— {l — r\)P where the choice 



TABLE 1. DSE results [10] for pseudoscalar and vector meson masses and electroweak decay 
constants, together with experimental data [11]. Units are GeV except where indicated. Quantities 
marked by t are fitted with the indicated current quark masses and the infrared strength parameter 
of the ladder-rainbow kernel. 





'"^1 = 1 GeV 


'"^ = lGeV 


- i^l)l=lGeV 


expt 
calc 


3 - 6 MeV 
5.5 MeV 


80 -130 MeV 
125 MeV 


(0.24 GeV)3 
(0.241 GeV)31' 
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expt 
calc 


0.138 0.131 

o.i38f o.or 


0.496 0.160 
0.497^" 0.155 


0.770 0.216 
0.742 0.207 


0.892 0.225 
0.936 0.241 


1.020 0.236 
1.072 0.259 



of rj is equivalent to a definition of relative momentum q; observables should not depend 
on rj. The meson momentum satisfies P^ = —M^. 

A viable truncation of the infinite set of DSEs should respect relevant (global) 
symmetries of QCD such as chiral symmetry, Lorentz invariance, and renormal- 
ization group invariance. For electromagnetic interactions and Goldstone bosons 
we also need to respect color singlet vector and axial vector current conserva- 
tion. The rainbow-ladder (LR) truncation achieves these ends by the replacement 

K{p,q;P) -)> -4;raeff(^^)^Jlv'^(^)|-7Ai ® t7v along with the replacement of the DSE 

kernel by Zig^D^y{k)V^ {q, p) -)> An tteff (^^) D^^f{k) fv y where k = p-q, and tteff (^^) 
is an effective running coupling. This truncation is the first term in a systematic ex- 
pansion [5, 6] of K; asymptotically, it reduces to leading-order perturbation theory. 
These two truncations are mutually consistent: together they produce color singlet 
vector and axial-vector vertices satisfying their respective Ward identities. This ensures 
that the chiral limit ground state pseudoscalar bound states are the massless Goldstone 
bosons from chiral symmetry breaking [7, 8]; and ensures electromagnetic current 
conservation [9]. 

We employ the ladder-rainbow kernel found to be successful in earlier work for light 
quarks [8, 10]. It can be written aeff(^^) = oc^^{k^) + a^^(A:^). The IR term implements 
the strong infrared enhancement in the region < k^ < I GeV^ required for sufficient 
dynamical chiral symmetry breaking. The UV term preserves the one-loop renormaliza- 
tion group behavior of QCD: aeff(^^) —^ as{k^y^°°'P{k^) in the ultraviolet with Nf = 4 
and Aqcd =0.234 GeV. The strength of Of^^ along with two quark masses are fitted 
to (qq), mj^i]^. Selected light quark meson results are displayed in Table 1. The in- 
frared kernel component is phenomenological because QCD is unsolved in such a non- 
perturbative domain. To help replace such phenomenology by specific mechanisms, it is 
necessary to first characterize its performance in new domains. 

Quark helicity and chirality in QCD are increasingly good quantum numbers at 
short distances or at momentum scales significantly larger than any mass scale. One 
manifestation of this is that, for chiral quarks, the correlator of a pair of vector currents 
is identical to the corresponding correlator of a pair of axial vector currents to all 
finite orders of pQCD. Non-perturbatively, the difference of such correlators measures 
chirality flips, and the leading non-zero ultraviolet contribution identifies the leading 
non-perturbative phenomenon in QCD. This is the four quark condensate [12]. Here, the 



vector correlator is formulated as the loop integral 

r-A 

'q 



n^v(^) = Jd'xe'''-'{0\Tj^{x)j+{Om = - j'^Tr{y^S{q+)Yl{q,P)S{q^)} , (3) 



where A indicates regularization, e.g., by the Pauli-Villars method, and Fy is the dressed 
vector vertex. The axial vector correlator is formulated in an analogous way and we di- 
rectly calculate the difference correlator which does not require ultraviolet regulariza- 
tion. We consider the transverse difference Yl\,^^{P^) = U^{P^) - nj.{P^). 

The leading non-perturbative contribution to D^ "^ starts with dimension d = 6 and 
involves the four-quark condensate in the form [13, 14] 



32nas{qqqq) .. a^iP^) 241 ..M^mi ^/ ^ 



Our numerical calculation of i'^n^^'^(i'^) identifies a leading ultraviolet constant rea- 
sonably well. The four quark condensate (qqqq) extracted via Eq. (4) is 65% greater 
than the common vacuum saturation assumption (qq)^ at the renormalization scale 
/i = 19 GeV used in this work. The low P^ limit provides a reasonable account of the 
first Weinberg sum rule [15, 16]: P^n^-^{P^)\p2^o = -fl in the fj, = 0.0924 GeV 
convention. This limit is due to Hy only and we obtain fji = 0.09 GeV that way. Our 
results are consistent with the second Weinberg sum rule [15] P'^Tlj^^{P^)\p2^^ = 0. 
The Das-Guralnik-Mathur-Low- Young sum rule [17] relates Jq dP^ P^U^^^{P^) to the 
strong component of mj^± — nij^n. We obtain 4.86 MeV for this mass difference in com- 
parison with 4.43 ± 0.03 from experiment. 

Data for the momentum-fraction probability distributions of quarks and gluons in 
the pion have primarily been inferred from Drell-Yan processes in pion-nucleon col- 
lisions [18, 19, 20]. For a recent review of nucleon and pion parton distributions see 
Ref. [21]. Lattice-QCD is restricted to low moments of the distributions, not the dis- 
tributions themselves [27]. Model calculations of deep inelastic scattering (DIS) parton 
distribution functions are challenging because it is necessary to have hard scale perturba- 
tive QCD features coexisting with a covariant nonperturbative model of soft scale bound 
states. In DIS calculations within the Nambu-Jona-Lasinio model [28, 29, 30], diffi- 
culties include a point structure for the pion and a marked sensitivity to the regulariza- 
tion procedure. In constituent quark models [31, 32], and instanton-liquid models [33], 
it is difficult to have pQCD elements join smoothly with nonperturbative aspects. An 
approach styled after the successful DSE treatment of the pion electromagnetic form 
factor [34, 35] can overcome these difficulties. 

The Bjorken kinematic limit of DIS selects the most singular behavior of a correlator 
of quark fields of the target with light-like separation z^ ~ 0. With incident photon 
momentum along the negative 3-axis, the kinematics selects z+ ~ z^ ~ leaving z^ as 
the finite distance conjugate to quark momentum component xP^, where x = Q^/2P ■ q 
is the Bjorken variable, q^ = —Q^ is the spacelike virtuality of the photon, and P is the 
target momentum. To leading order in the operator product expansion, the associated 
probability amplitude qf{x), characteristic of the target, is given by the correlator [36, 
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FIGURE 1. Left Panel: Pion valence quark distribution evolved to (5.2 GeV)^. Solid line is the full DSE- 
BSA calculation [22]; dot-dashed line is the semi-phenomenological DSE-based calculation of Hecht et 
al. [23]; experimental data points are from [20] at scale (4.05 GeV)^; the dashed Une is the recent NLO 
re-analysis of the experimental data [24]. Right Panel: The ratio of u-quark distributions in the kaon and 
pion. The solid line is our preliminary result from DSE-BSE calculations [25, 22, 21]; the experimental 
data is from [26, 18]. 
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(5) 



where / is a flavor label. Here the 4-vector components that arise naturally are 
a"^ — {a'^ ±a^)/\/2. This probability amplitude is invariant, it can easily be given a 
manifestly covariant formulation, and its interpretation is perhaps simplest in the infi- 
nite momentum frame where qf{x) is the quantum mechanical probability that a single 
parton has momentum fraction x [38]. Note that qf{x) = —qA—x), and that the valence 

quark amplitude is q'^Ax) = qf{x) —qr[x). It follows from Eq. (5) that /q dxq^'Ax) = 
{7r{P)\jf {0)\7r{P)) /2P^ = Fjj:{0) = 1. Approximate treatments should at least preserve 
vector current conservation to automatically obtain the correct normalization for valence 
quark number. 

In a momentum representation, qf{x) can be written as the evaluation of a special 

Feynman diagram [36, 39] qf{x) = y^d{k+ -xP+) ticdir^ G{k,P)] where G{k,P) 
represents the forward ^-target scattering amplitude. Ladder-rainbow truncation, which 
selects the valence qq structure of the pion, yields 



^/(. 



X 



tv,d[r^{p,p)s{p)r+ip-x)s{p)r„{p,p)s{p-p)] 



(6) 



where trcd denotes a color and Dirac trace, and r+(p; jc) is a generalization of the dressed 
vertex that a zero momentum photon has with a quark. It satisfies the usual inhomoge- 
neous BSE integral equation (here with a LR kernel) except that the inhomogeneous 
term is 7+ d (/>+ — xP^). This selection of LR dynamics exactly parallels the symmetry- 
preserving dynamics of the corresponding treatment of the pion charge form factor at 



q^ = wherein the vector current is conserved by use of ladder dynamics at all three 
vertices and rainbow dynamics for all 3 quark propagators [34, 35]. Here the num- 
ber of valence M-quarks (or d) in the pion is automatically unity since the structure of 
Eq. (6), along with the canonical normalization of the qq BS amplitude rjt{p,P), ensures 
/q dxq^{x) = 1 because /q dxr^{p;x) gives the Ward Identity vertex. 

Eq. (6) is in Minkowski metric so as to satisfy the constraint on p+, but LR dynamical 
information on the various non-perturbative elements such as S{p) and Yji{p,P) is 
available only in Euclidean metric [10]. Since qf{x) is obtained from the hadron tensor 
W^^ which in turn can be formulated from the discontinuity T^^{£) — r^^(— e), we 
observe that all enclosed singularities from the difference of Wick rotations cancel 
except for the cut that defines the object of interest. With use of numerical solutions 
for dressed propagators and BS amplitudes, that give an accurate account of light 
quark hadrons, our DIS calculations significantly extend the exploratory study made in 
Ref. [23]. That work employed phenomenological parameterizations of these elements. 

In Fig. 1 we display our DSE result for the valence w-quark distribution evolved to 
2^ = (5.2 GeV)^ in comparison with kN Drell-Yan data [20] with a scale quoted as 
2^ > (4.05 GeV)^. We also compare with a recent NLO reanalysis of the data at scale 
Q^ = (5.2 GeV)^. The distribution at the model scale Q^ is evolved higher by leading 
order DGLAR The model scale is found to be Qq = 0.57 GeV by matching the x'^ 
moments for n = 1, 2, 3 to the experimental values given independently at (2 GeV)^ [40]. 
Our momentum sum rule result /g dxx{u-n; + d-n) = 0.74 at Qq clearly show that in a 
covariant approach the retardation effects of one gluon exchange assign some of the 
momentum to gluons. The corresponding momentum sum for the kaon is 0.76. 

The ratio uk/uti measures the dynamical effect of the local environment. In the kaon, 
the M-quark is partnered with a significantly heavier partner than in the pion and this 
shifts the probability to relatively lower x in the kaon. Our preliminary DSE model cal- 
culation [25, 22, 21] is shown in Fig. 1 along with available Drell Yan data [26, 18]. 
Here we include only the leading two invariants of the pion BS amplitude, E{q^P) 
and F{q,P), where q is qq relative momentum. For both amplitudes only the lowest 
Chebychev moment in q- P is employed . This variable does not occur in static quan- 
tum mechanics, nor in the Nambu-Jona-Lasinio point-coupling field theory model [28] 
which also neglects the q^ dependence. We do not make such a point meson approxima- 
tion here; the q^ dependence comes from the BSE solutions. Nevertheless, the essential 
features of the ratio uk/u-h are adequately reproduced by a generalized Nambu-Jona- 
Lasinio model [21]. 
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